#include "cCoupledsq.h"

Coupledsq::Coupledsq()
{
	setName("Coupled finite sq. well");
		setUniverse(40.0);
		setWidth(6.0);
		//setPotLeft(-1.0*getWidth()/2.0);
                setPotLeft(-1.0*getUniverse()/2.0);
		//setPotRight(getWidth()/2.0);
                setPotRight(getUniverse()/2.0);
		setFx(getUniverse()/2.0);
		setTx(-1.0*getUniverse()/2.0);
                setUniLeft(-1.0*getUniverse()/2.0);
                setUniRight(getUniverse()/2.0);
		setDepth(5.0);
		setMinEnergy(-1.0*getDepth()+0.05);
		//filter[activefilter]->calculate(this);
                temp.push_back(-1.0 * getWidth()/2.0);
                temp.push_back(getWidth()/2.0);
                temp.push_back(-0.5);
                temp.push_back(0.5);
                setDiscontVector(temp);
}

double Coupledsq::getPotentialValueAt(double x) {
	if(fabs(x) > getWidth()/2.0) {return 0.0;}
	else if(fabs(x) <= 0.5) {return 0.0;}
	else {return fabs(getDepth())*-1.0;}
}

Coupledsq::~Coupledsq()
{
}

/*  Coupled finite square well
 * 
 * _____        ________	_______
 * 	|	|	|	|
 * 	|	|	|	|
 * 	|	|	|	|
 * 	|_______|	|_______|
 */ // -3   -0.5       0.5      3
